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Abstract 

The nonlinear oscillations of a scalar field are shown to have cosmological equations 
of state with w = p/p ranging from — 1 < w < 1. We investigate the possibility that 
the dark energy is due to such oscillations. 
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Astrophysical data support the existence of dark energy [1, 2}. Since many proposed 
solutions of the cosmological constant problem lead to exactly zero vacuum energy for empty 
space, it is natural to consider so-called quintessence models in which the dark energy is 
comprised of some scalar field which is slowly evolving towards its minimum [2]. The main 
objections to these models are their typically unnatural potentials, and that they require the 
suppression of higher dimension operators likely to be induced by quantum gravity [3] . 

In this letter we investigate a qualitatively different idea: that the dark energy is due to 
(possibly rapid) nonlinear oscillations rather than slow evolution on cosmological timescales. 
We consider oscillations in scalar models 

C=±(d>z) 2 -V(z), (1) 

where the potential V(z) = a\z\ l near its minimum. Potentials with / < 2 are particularly 
interesting, as we will see below that they yield an equation of state w = p/p < 0. This form 
of V(z) may appear odd, but a change of variables to, e.g., = \z\ 1 / 2 yields 

C= l -K{<t>){d^) 2 -a<t>\ (2) 

with K((f>) = (j^j ifi^r 1 . A kinetic term of this type can be obtained from the Kahler 
potential in supersymmetric models. We focus here on the classical behavior of z, but its 
quantization when I < 2 warrants further investigation. It has the unusual property that 
there are no perturbative degrees of freedom - that is, small oscillations about z = have 
infinite frequency, since V"(z = 0) does not exist. Only large (non-infinitesimal) z oscillations 
can have finite energy density. This may lead to a number of interesting features: one might 
expect that z decay and production rates, as well as radiative corrections, only arise from 
nonperturbative effects and are exponentially small. 

The redshift of a field undergoing nonlinear oscillations can be calculated through its 
average equation of state and depends on the ratio w — p/p. From the scalar field equations 
in a Robertson- Walker universe, one obtains 

/ D \ 3(1+10) 

where R is the Robert son- Walker scale factor. The scalar field energy redshifts like radiation 
when w — 1/3, like matter when w — 0, etc. 

It remains to calculate the equation of state for nonlinear oscillations. We note that 
p = T — V and p = T + V, where T is the kinetic energy density and V is the potential 
energy density. We calculate the relation between T and V averaged over an oscillation 



2 



period, which is much smaller than the cosmological timescale except during the very first 
oscillations, which begin when the age of the universe is of order V"(z) . 



We define 



(T) = \ f dt i 2 (4) 



and 



(V) = Jdtz 1 , (5) 

where each integral is taken over the same period with boundary conditions i(0) = z(t) = 
(or equivalently z = at the endpoints) and we have adopted units in which the overall scale 
of the potential is unity. Since 

— (zz) =z 2 + zz , (6) 
and the equation of motion is z — —Iz 1-1 , we can rewrite the average potential energy as 

(V)= 1 1 Jdtz 2 = j(T) . (7) 



This yields 



w 



(1-2) 



(1 + 2) ' 

with — 1 < w < 1. In potentials with I < 2, the average potential energy dominates over the 
kinetic part, and the pressure is negative. In the limit I — > these oscillations behave like a 
cosmological constant. In higher order potentials (/ > 2), the situation is reversed, leading 
asymptotically to w = 1 as / — > oo , or p ~ 1/i? 6 . These oscillations redshift away rapidly, 
although it was noted in [4] that the large 1 behavior of w — 1 can never be achieved, due 
to an instability to a non-oscillatory scaling solution. 

Given a periodic solution to the z equations of motion, one can obtain a rescaled solution 

via 

z(t) -> z (at) . (9) 



For I < 2, the frequency of oscillation goes to infinity as the amplitude goes to zero. Note, 
however, that the average energy density goes to zero in this limit. 

An advantage of nonlinear oscillation models of quintessence is that the potential V(z) 
need not be characterized by the size of the current dark energy density ~ (10~ 3 eV) 4 , 
nor be fine-tuned to be flat (i.e., have curvature of order the inverse horizon size squared 
~ (10 -33 eV) 2 ). Rather, the potential can be characterized by a larger energy scale more 
familiar to particle physics, with no small dimensionless parameters. The smallness of the 
energy density today relative to this scale could be explained by a small z oscillation ampli- 
tude. 
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One scenario that would result in a small oscillation amplitude is if the original energy 
density in the z field were diluted away by inflation, and the subsequent reheat temperature 
insufficient to repopulate it. This is quite plausible if the couplings between the inflaton and 
ordinary matter fields to the z field are small, for example suppressed by the Planck scale 
if z is a hidden sector field 1 . The current z energy density is dependent on the number of 
e-foldings N during inflation 



where v = 3(1 + w), T m< j ~ 5 eV is the temperature at which the universe becomes matter 
dominated and T r h is the reheat temperature after inflation. For example, using v = 1/2 (or 
w = -5/6, consistent with the WMAP limit of w < -0.78 [1]), T rh = 5- 10 10 GeV and p, t the 
Planck energy density, we find that N ~ 510 in order that ptoday ~ (10 _3 eV) 4 . For smaller 
Pi ~ (10 11 GeV) 4 , appropriate for intermediate scale inflation, we find iV ~ 370. Note, while 
the scenario described here explains the small dark energy density today, it does not address 
the question of coincidence: why is p z of order p C riticai today? 

If the energy scale characterizing V(z) is small (within several orders of magnitude of an 
electron volt) little dilution is necessary, and it could be provided by the expansion of the 
universe after inflation. If z originates from the scalar component of a chiral superfield $, 
the energy scale of V(z) (i.e., the parameter a in (2)) is protected by supersymmetry from 
radiative corrections. A Yukawa coupling between z and its superpartner can be excluded 
by imposing $ — > — $ symmetry in the superpotential, thereby stabilizing z. 

To conclude, we find that realistic quintessence models based on nonlinear oscillations 
can be constructed without any fine-tuning of fundamental parameters. These models will 
be disfavored if future data show that w — — 1. 

Note added: After this work was completed, we learned that the result (8) for the equation 
of state was previously derived by M. Turner in [5]. Also, models based on a hyperbolic cosine 
potential raised to a fractional power have been considered in [6] , which at late times exhibit 
oscillations of the type considered here. 



^ome reheating of the z field is inevitable, even if its couplings to the inflaton are very small. However, 
the thermal z bosons produced do not necessarily contribute to the coherent oscillations studied here - their 
energy rcdshifts away more rapidly. 




(10) 
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